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We evaluate the operator product expansion (OPE) for a mixed correlator of the isovector and 
isoscalar vector currents in the background of the nucleon density with intrinsic isospin asymmetry 
[i.e. excess of neutrons over protons] and match it with its imaginary part, given by resonances and 
continuum, via the dispersion relation. The leading density-dependent contribution to p — uj mixing 
is due the scattering term, which turns out to be larger than any density dependent piece in the OPE. 
We estimate that the asymmetric density of n n —n p ~ 2.5 x 10 -2 fm 3 induces the amplitude of p — u) 
mixing, equal in magnitude to the mixing amplitude in vacuum, with the constructive interference 
for positive and destructive for negative values of n n — n p . We revisit sum rules for vector meson 
masses at finite nucleon density to point out the numerical importance of the screening term in the 
isoscalar channel, which turns out to be one order of magnitude larger than any density-dependent 
condensates over the Borel window. This changes the conclusions about the density dependence of 
m u , indicating ~ 40 MeV increase at nuclear saturation density. 



I. INTRODUCTION 



Changes of hadronic properties in hot and dense nuclear medium are an intriguing issue which ties together modern 
particle and nuclear physics. The interest to these questions has been intensified over the past decade due to the 
possibility of studying the transition from hadrons to the deconfining phase at heavy ion collisions. In particular, the 
modification of vector meson properties in nuclear medium has been a subject of a persistent theoretical activity [EJ. 
This was initiated by the idea that in nuclear medium the vector meson masses should drop as a precursor to the 
chiral symmetry restoration pi. Several experiments have also been proposed to study the changes of masses, widths 
and coupling constants of vector resonances in dense (and/or hot) nuclear matter ||. 

The properties of vector resonances in vacuum and the effects of isospin symmetry violation on the mixing of the p, 
uj resonances in vacuum have been investigated rather carefully by means of QCD sum rules in the past In the 

pioneering work of Ref. it was found that the nonzero value for the p — uj mixing can be linked to the difference 
of light quark masses, and the possibility of m u — is seemingly excluded. 

Later, the QCD sum rule method was extended to finite temperatures and densities jjj. A number of analyses 
|p|--|lT| have found that the masses of p and uj resonances decrease in nuclear medium]^]. In Refs. @,^5) finite widths 
of the vector mesons have been taken into account by hand and by calculation of the r/io-meson self-energy in a 
chiral model for the spectral function, respectively. For the p-meson channel it was found in Ref. [|l4| that at nuclear 
saturation density an increasing width of the /^-resonance necessitates an increasing p-meson mass. However, for large 
values of the width the mass is blurred over a large window of possible values. 

While appreciable efforts have been directed to estimate the density dependent modification of the masses and 
lifetime of the light vector mesons at finite density (and/or temperature), the question of p — uj mixing at finite 
densities (and/or temperature) has not received much attention. In fact finite nuclear densities can have a significant 
impact on this amplitude. The fact that nuclear matter can intrinsically be isospin asymmetric implies that the p — uj 
mixing in matter can potentially be larger than the vacuum part of the mixing which is induced by the difference in u 
and d quark masses, small in units of characteristic hadronic scales. This idea was suggested first in Ref. |l6) where 
it has been pointed out that the presence of asymmetric nuclear matter has a profound effect on the mixing of the p 
and uj resonances. There the mixing angle was determined from the matter induced non-diagonal self energy of the 
p resonance by employing an SU(2) p symmetric hadronic model. Subsequently such a matter induced mixing has 



1 See, however, the work of Y. Koike 113] , where opposite behaviour is claimed. A later analysis, based on the relation between 
the current-nucleon forward scattering amplitude and the scattering length of the vector meson off the nucleon in the static 
limit, again revealed negative mass shifts in the linear density approximation |13[. 
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also been analysed on a more elaborate footing in Ref . |L7| . Along the same lines the author of Ref. [H investigated 
the nucleonic density and temperature dependent p°-oj mixing at a fixed asymmetry. Thereby, an enhancement of the 
modulus of the vacuum mixing amplitude was found due to finite density. In all model descriptions the vacuum part 
of the mixing serves as an input parameter, to which all the results are normalized by hand in the limit of vanishing 
density. To this end, it is desirable to obtain an independent analysis of the mixing using finite density QCD sum 
rules, which hopefully would allow to treat the vacuum part and the density part from the first principles. At this 
point we already notice one principal problem of finite density QCD sum rules. In the presence of nuclear matter 
there exist nonscalar condensates which can be related to the twists of different dimension. In general, going from 
mass dimension 2n to 2(n + 1) the ratio of contributions with a nonzero, fixed twist 2k is 



This requires the external momenta to be much larger than for the OPE to converge. However, the possibility to 
link properties of a ground state resonance to nonperturbative effects in the vacuum (the condensates) via the sum 
rule requires external momenta of ~ m^r. We will later show that for the contribution of twist operators there is a 
numerical suppression in the corresponding Wilson coeffcients up to mass dimension six. Since at higher dimensions 
we have no parametrical smallness the above should limit the applicability of QCD sum rules at finite nucleonic 
density. 

In this paper we study the behavior of the isoscalar-isovector mixed correlator of the two vector currents in order 
to extract nuclear density effects. The asymptotic behaviour of this correlator at large space-like external momenta 
can be studied within the perturbative QCD framework, with the power corrections represented by quark, gluon, 
quark-gluon, etc. condensates. In the presence of finite nucleon density the power correction due to these condensates 
will change as compared to their vacuum values. Due to the presence of the preferred reference frame, in which the 
nucleons are at rest, new density-dependent power corrections will appear. In both cases we assume the small density 
regime and keep only the linear terms in the external nuclear density. As we shall see, this approximation is justified 
for densities not larger than the nuclear saturation density, which is small in proper "vacuum" units. 

The asymptotics of the two-point correlation function, calculated this way, can be related to the "phenomenological 
part" which includes the contributions of vector resonances, continuum and the screening terms A success or a 
failure of the QCD sum rule analysis of vector meson properties would depend on how reliably the contribution of 
individual resonances (p, a;,...) can be separated from the rest of the contributions. 

We carefully examine the density-dependent part of the operator product expansion (OPE) and find that the effects 
of matter-induced mixing due to nucleonic matrix elements of nonscalar and scalar QCD operators in asymmetric 
nuclear matter follow a certain hierarchy. The asymmetric density-induced effects start dominating vacuum contri- 
bution at asymmetries a pn = (n p — n n )/(n p + n n ) w 0.2 and an overall nucleonic density twice the nuclear saturation 
density n° N = 0.17 fm -3 =(111 MeV) 3 . However, the analysis of the phenomenological part of the QCD sum rules 
shows that the scattering contribution, usually called screening term, turns out to be numerically by far more impor- 
tant. Brought to the OPE-side of the sum rule, the screening term can be regarded as a mass dimension two power 
correction. Already at intermediate asymmetries a pn ~ 0.1 and saturation density vacuum and matter induced p-w 
mixing are of the same sign and comparable in magnitude. 

The smallness of the density-dependent pieces in the OPE as compared to the screening term indicates that any 
conclusion about the density-dependent piece in the p — u> mixing amplitude will mostly depend on the assumptions 
made about the spectral density, i.e. what is usually called phenomenological part of the sum rules. This casts strong 
doubts on the applicability of the finite density QCD sum rules for the extraction of the isovector-isoscalar mixing 
since the density-dependent "QCD input" is negligibly small. This concern lead us to re-examine the screening terms 
in the isovector-isovector and isoscalar-isoscalar correlators which were used in previous works to investigate 

the modification of p and lo masses in nuclear matter. We have found that all previous analyses have used the same 
value for the screening terms in the isovector-isovector and isoscalar-isoscalar correlators. This is an unfortunate error 
because the screening term in the omega channel turns out to be 9 times larger than the value used in Refs. p|-pT|. 
This changes drammatically all the conclusions about the behaviour of m u in nuclear matter, and indicates that m w 
is a growing function of density in the linear density approximation. 

II. ISOSINGLET-ISOTRIPLET CORRELATOR OF THE TWO VECTOR CURRENTS AT FINITE 

DENSITIES 

We start with the (causal) mixed correlator of isotriplet and isosinglet currents in asymmetric nuclear matter 
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where 



= d*xe^(TjT(x)jS(0)) nN , 



J'J = \ ( u Jv u ~ d ln d ) > if = \ ( u Jv u + d lv d ) 



(2) 
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We choose the same normalization of the two currents, which also means that their couplings to physical p and uj 
resonances are approximately equal. In Eq. (||) the Gibbs average ( } njv (tin indicating finite nucleon density) is 
approximated by a vacuum and one-particle nucleon states p|-|ll|] . Due to the presence of a singled out rest frame 
with four-velocity u M there are, in general, two independent, current conserving tensor structures (longitudinal and 
isotropic) into which can be decomposed. However, in the limit <f — > one of the corresponding invariants II;, H 
becomes redundant B, and we therefore concentrate on 11/ which satisfies the following dispersion relation [p| 
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Subtracting the terms attributed to the p — > 7 — ► ui electromagnetic mixing from the spectral representation B and 
the OPE and appealing to the literature on density dependent OPE's of p and uj current correlators we arrive at the 
following sum rule 
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The asymptotic behaviour of the lhs of Eq. (g) can be calculated by means of the operator product expansion (OPE) . 
The result is given in terms of the perturbative contribution and power corrections, proportional to the condensates, 
taken in the presence of the external nucleon density. Retaining terms up to the order Q^ 6 , we present the result in 
the following form: 
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In Eq 




h, symbol S denotes the operation of making tensors symmetric and traceless. As usual, (for example 
| Jl2]Jl5[ |), the averages over mixed operators and twist four contributions have been omitted in Eq. (^|). The 
former can either be reduced to four quark operators by use of the equation of motion (these contributions are already 
included), or they are suppressed at p 2 « 1 GeV 2 since there the gluon content of the nucleonic wave function is 
small ||. The latter has been argued in Ref. || to have no substantial effect on the p and ui mass shifts, and we 
will therefore omit twist four operators. Further progress in calculating the OPE depends on how accurately we can 
predict the size of various contributions to II'. We restrict ourselves to the case of low and medium densities, so 
that the linear (mean field) approximation is justified, and the density dependent part enters in the final expression 
multiplied by the matrix elements over the single nucleon states. We further make use of the vacuum saturation 
hypothesis which becomes an exact relation in the limit of large number of colors. This hypothesis is known to 
"work" reasonably well in vacuum. However, its application to the nucleon matrix elements is not fully justified. 
We use this hypothesis to estimate the order of magnitude of dim 6 contribution, noting that their numerical weight 
in the final result turns out to be small as the OPE is largely dominated by dim 4 contributions. With all these 
assumptions, Eq. (o) can be reduced to the following form: 
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In Eq. (Q) njv (%>, denotes the total nucleonic (proton, neutron) density; a pn the p-n asymmetry, defined as 
a pn = {n p — n„)/njv; m p the proton mass; {qq} the value of the u-quark condensate; 7 the asymmetry of u-quark and 
ci-quark condensate defined as 7 = (dct) / (uu) — 1; E^jy = (45 ±7) MeV || the nucleon sigma term; and to is defined 
as to = l/2(m u + m ( i). The electromagnetic coupling is a, the strong coupling at scale fi is a s (fi 2 ). In Eq. (0) we have 
already used the numerical smallness of isospin and chiral symmetry violating parameters as compared to the normal 
hadronic scale and thus neglected terms proportional to m„(d)7, 7 2 and so on. For similiar reasons it is justified to 
neglect the effects of isopsin violation in nucleon matrix elements and take (p\uu — dd\p) = — (n\uu — dd\n), which 
leads to the dependence on the asymmetry factor a pn . The scalar matrix element (p\uu — dd\p) is related to the 
baryon octet mass splitting, (to- — m^)/m s , and it is numerically close to 0.7 |l£|. 

As for the contribution of symmetric and traceless twist two quark bilinears [ 3[ of dimension four and six, their 
nucleonic matrix elements are determined by the quark parton distributions A\~ (/x 2 ) and A^~ d (ii 2 ) as || 
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where k^ denotes the nucleon momentum, D p is the gauge covariant derivative. In general, the factor A q k ([i 2 ) can be 
obtained from the parton distributions Q(x,fi 2 ) and Q(x,fi 2 ) in the proton as 
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In Ref. pl| the parton distributions in the nucleon have been fitted to experiment at a resolution scale /i 2 
GeV 2 . Using these distributions and performing the integrations of Eq. (|J), we obtain 



A%- d {fi 2 = 0.26 , GeV^) = 0.429 , A\~ a {^ = 0.26 , GeV^) = 0.097 . 
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To generate the respective values at the scale /z 2 « 1 GeV 2 relevant for the sum rule we simply use the conversion 
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at 1 GeV from Ref. 



we arrive at 



^2 _d (l GeV 2 ) = 0.32 , A u - d (l GeV 2 ) = 0.062 . 
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Keeping this in mind, we will for now proceed to the numerical evaluation of OPE. Performing a Borel transformation 
of the OPE of Eq. (Q) and omitting the numerically strongly suppressed dimension two power correction, we obtain 
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In this expression we have used the following set of values: 

m u = 5MeV; m d = 9MeV; to = 7MeV; m p = 940 MeV; n N = n° N = (HlMeV) 3 (the nuclear matter saturation 
density); E„ N = 45 MeV; {qq} = -(225MeV) 3 ; 7 = -10~ 2 @; a = 1/137; and a s (lGeV 2 ) = 0.5 §. The quark 
condensate has been factored out numerically for the sigma-term and the twist contributions. 

Several important observations should be made at this point. For M ~ 1 GeV, a pure perturbative contribution is 
negligibly small as compared to power corrections. The latter are dominated by dimension 4, with the constant term 
originating from m u and md mass difference and the a np njv-dependent piece given by the A^ - ^ contribution. At the 
level of dimension 6 we observe three different terms (second line of Eq. (O)): vacuum part, density-dependent scalar 
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condensate and twist contributions. At this dimension the density dependent contribution from scalar condensate and 
twist tend to cancel each other. This cancellation can be an artefact of chosen parameters and/or of the crude nature 
of approximations made in estimating the size of the four-quark matrix elements over the nucleon. Nevertheless, 
M ~ 1 GeV, and the OPE is dominated by dim 4 terms, where at ^ff ~ 1 the suppression of dim 6 is about 50%. 

For higher values of asymmetric density II' (M 2 — 1 GeV) changes sign. 

What does this behaviour of II' (M 2 = 1 GeV) mean in terms of the u> — p resonance mixing amplitude? To answer 
this question we should parametrize the spectral function in terms of the resonance contributions and analize the 
resulting sum rule (|). 

Following Refs. pJ7|, we approximate the imaginary part of the correlator by contributions of p, u>, p' , ui' resonances 
and continuum: 
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The contribution to the mixing due to the electromagnetic continuum is small Q. Therefore, we will neglect it in the 
subsequent consideration. In Eq. ( |l3|) f p and f u refer to the p and to residues of the p-oj current propagator, and p' 
and uj' symbolize the cumulative effect of higher resonances^ introduced in the original analysis || in order to have 
consistent asymptotic behaviour of II' (M 2 ). Besides the "usual" annihilation continuum above a certain threshold soj 
Eq. ( |l3| ) exhibits a scattering term which behaves as a pole at s = (Landau pole) |?J . The corresponding coefficient 
can be calculated explicitly, and the result in the leading order in Fermi momentum (pf/rriN expansion) is given by 

P S sI = — [FpF^n p + F^F^n n ] = —a pn n N . (14) 

Here the coefficients ar ^ denned via nucleon matrix elements of quark vector currents at vanishing momentum 

transfer: 

Fp = F* = (p\u lQ u + d l0 d\p) = 3 
Fp = -F^ = (p\u l0 u - d l0 d\p) = 1. (15) 

After Borel transformation the contribution of the Landau screening term is usually carried to the lhs of the sum rule 
to effectively become a power correction of dimension two in the expansion of IIJ(M 2 ). Defining f pu = 1/2 (f p + f u ), 
to 2 = 1/2 (m 2 + m 2 ), Am 2 = m 2 — m 2 and fi = (f u — f p )fh 2 /(f pLJ Am 2 ) (the primed quantities are defined 
analogously), we quote the result of Ref. B relating f puJ to the measurable quantities to 2 , Am 2 , g p , and g u 
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where g p , g u are the respective decay constants, and 8 pbJ enters the measurable mixing parameter e as follows 
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Thereby, e is defined as 

to = uq + spa , p = po - eu>Q , (18) 

and r denotes the width of the respective resonance. It is fair to remark at this point that the observable combination, 
e ~ Sp^T^rrip 1 will have an additional dependence on density due to a substantial increase of T p with tin p3[ . Thus, 
finding the decrease of £ with density would certainly allow to conclude that e is decreasing. The opposite behavior, 
a rising £, would complicate the prediction of e(ri7v). 



In Ref. |4| the cumulative values m p /,m^; were chosen to be about 1.5 GeV , which is well below the physical masses (~ 1.7 
GeV) of the resonances p',uj' . 
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The final sum rule is given by the following expression: 



" . ,0-G.V- $ + £ ( 2 MeV- ^MeV H ^ ] - 



where all masses and the Borel parameter are taken in units of GeV. It is remarkable that the screening term, brought 
to the OPE side of this sum rules, completely dominates other density dependent contributions. This shows that in 
the asymmetric nuclear matter background the influence of the screening term on the u — p mixing is by far more 
important than any changes of the QCD condensates. Moreover, for any realistic M 2 the screening term becomes 
comparable to the vacuum contribution to the mixing at tin — n N and asymmetries as low as a np ~ 0.05. 

In the limit of vanishing density, relation ( [l9|) reduces to the known sum rule for p — ui mixing. A naive evaluation 
of this sum rule at p-meson mass, M 2 = (0.77) 2 , and at tin — 0, gives a reasonable agreement with experimentally 
measured value £ = 1.1 X 10~ 3 p2) with (3 ~ 0.5, advocated in Ref. ||. Next, we parametrize the linear dependence 
of £ and (3 on the density as follows: 

£ = e (o)+£Wn; = 0W+0Wn, (20) 

where n denotes a np riN in units of 0.2n^. 

The primary reason for the introduction of the p' — u>' contribution in Ref. |6) was the absence of the 1/M 2 term 
in the OPE side of the sum rule, so that p and u> contribution alone would not be consistent with the asymptotic 
behaviour of II'. Thus, the role of p' — w' is to imitate the cancellation of 1/M 2 terms in contributions of various 
resonances at large M 2 . For a semiquantitative determination of the linear density dependence of £ and (3 we proceed 
as in Ref. Q . There the vacuum values of £ and (3 were estimated by choosing M = m p , which strongly suppresses 
the higher resonances. It should then be legitimate to compare powers of M~ 2 in the OPE and the lowest resonance 
contribution. The result is given by 

2 l Z — -2 0-10 -4 - (2\) 

£(0) + 0(0) ^ AU £(03/3(0)^2 • ^ Zl > 

Using this relation, we can find (3^ and £W separately, evaluating ( |l9| ) at M 2 — 0.59. The final estimate of 
reads as 

~ [2.3 - 0.8] x 10" 3 = 1.5 x 10" 3 , (22) 

where 2.3 originates from the screening term and -0.8 comes from the OPE. A similar number can be obtained from 
the combination of ( |l9| ) and its first derivative in M 2 . This value of ^ < - 1 - 1 leads to the doubling of mixing amplitude 
and complete screening at n n — n p ~ ±2.5 x 10~ 2 fm 3 , respectively. 



III. IMPORTANCE OF THE SCREENING TERM FOR THE ISOSCALAR-ISOSCALAR CORRELATOR 

Having found such an important role of the screening term in the isoscalar-isovector mixed correlator, we would 
like to return to previous analyses of diagonal correlators (isovector-isovector and isoscalar isoscalar) which were used 
to extract the behaviour of m p and at finite nucleon density In all these papers it was found that masses 

and coupling constants of p and to resonances behave similarly in nuclear matter, simply because the OPE sides of 
the sum rules in both cases are the same after the application of the vacuum saturation hypotheses. 

We use the same symmetric normalization of the two currents, Eq. ([}]). From now on we neglect the asymmetry of 
the nuclear matter and other isospin breaking effects. Then the sum rules for isovector-isovector and isoscalar-isoscalar 
correlators in medium take the following symbolic form: 
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where C4 and cq are the same for both expressions. Obviously, at vanishing nuclcon density F u 



It is remarkable, that the screening terms in Eqs. 
screening term in the isoscalar-isoscalar channel is due to 



F p and S u 



are different by a factor of 9. The enhancement of the 



ss 

PSc 
TT 
PSc 



p v _ q 

pTpT 
P P 



(25) 



This difference was overlooked in Refs. |8|,[|,[n] [j]. 

The coefficients C2 and C4 can be computed along the same standard technique (again with considerable degree 
of uncertainty for cq). When plugging these values into the sum rules ( p3f|24| ), we obtain the following numerical 
relations: 
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where again all masses and dimensional coupling constants are taken in GeV units. It is remarkable that at M ~ lGeV 
the screening term in the to sum rule is larger by an order of magnitude than any other density-dependent term from 
the OPE! 

As in the previous case, it is convenient to parameterize the density dependence of masses and coupling constants 
as follows: 
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Using the sum rules (g6() and (^7|), an( ^ the first derivatives of these expressions, we solve for to* 2 as a function 
of Sq, and Borel parameter M. The dependence of the threshold on the the density is obtained by requiring the 
Borel curves, m* (M 2 , S* , n), be parallel over the Borel window which we take from 0.6 to 1.2 GeV for different 
values of densities. The slope of the Borel curve m(M 2 ) in the Borel window at zero density represents a "systematic 
uncertainty" introduced by sum rules and the requirement of the Borel curves to be parallel at different densities is 
equivalent to the requirement that this uncertainty does not change while going to finite but small densities. The 
resulting dependence of So on the density, Sq 1 ^ /Sq ^ = —0.2 for p and —0.1 for u), allows us to deduce the following 
estimates for the linear dependence of masses on the density: 
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Our estimate for m p agrees with the results of previous analyses [pp|Jll[|. The result for mi/ has the opposite 
sign and correspond to an 40 MeV increase of m w at the nuclear saturation density. This difference could be easily 
explained by the error in the screening term for u> sum rule in ^,^,|llj . The disagreement with the results of |l5| ] , where 
the correct form of the screening terms is used, is harder to explain, and we hypothesize that it could be an artefact 
of different numerical methods used to extract the density dependence of the resonance masses and thresholds. 



IV. DISCUSSION 



Apart from the question of (non)convergence of the OPE we would like to point out some concerns about usefulness 
and validity of the sum rules at finite densities. 



3 In Ref. |llj p ac was taken as a free search parameter and determined from the sum rules at the level consistent with pgj for 
both correlators. It casts a strong doubt on the validity of the whole approach, since the actual value of the screening term for 
to should be 9 times larger. 
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1. Vacuum factorization at dim 6. It is unclear what the status of factorization procedure is, especially in the 
presence of nuclear matter. In principle, one could try to relate four-fermion matrix elements over the nucleon 
states, which appear in the calculation of the OPE, to some measured processes induced by weak interactions. 
Indeed, non-leptonic hyperon decays and parity violating pion-nucleon coupling constants could be reduced to 
similar matrix elements from the four-quark operators. It is unclear, though, whether such an analysis is feasible. 

2. The importance of a particular choice of the spectral function. In linear density approximation the analysis of 
the examples of the p — u and u — u sum rules suggest that there are large contributions from the respective 
screening terms. In fact, these contributions dominate all density dependent pieces in the OPE. It means that 
the "QCD input" in these channels is not important in comparison with the choice of the spectral function at 
finite density. 

3. Is the linear density approximation valid up to n and beyond? The use of the dilute Fermi gas to model the 
behaviour of the scattering terms and the QCD condensates has its limitations, and a more realistic description 
may greatly affect the resulting sum rule. However, it seems unfeasible to calculate QCD operator averages 
over interacting multi nucleon states which one would have to consider when going beyond the dilute gas 
approximation. An inclusion of Fermi motion of noninteracting nucleons is practically doable but does not alter 
the zero momentum result significantly. 

In conclusion, we have considered the correlator of isovector and isosinglet vector currents in the presence of the 
asymmetric nuclear matter in linear density approximation. We see a significant dependence of the OPE on n n — n. pi 
which becomes comparable to vacuum contributions at n n — n p ~ 0.05n^. An attempt to extract the p-uj mixing, 
using the dispersion relation has shown that this mixing is more affected by the presence of the scattering term 
than by density-dependent part of the OPE. A similar tendency exists in the isosinglet-isosinglet channel, which is 
normally used to deduce the dependence of m w on density. Hence, in linear density approximation the explicitely 
density dependent part of the spectral functions (scattering terms) in the p-ui and w-lo channels dominantly drive the 
density dependence of hadronic parameters. The density dependence of the width of the resonances, which has been 
neglected here, does not alter this finding. However, it may drastically change the conclusions about the direction of 
resonance mixing at finite, asymmetric density. 



ACKNOWLEDGMENTS 



R.H. and M.P. would like to thank V. Eletsky and A. Vainshtein for sharing their scepticism about the relevance 
of QCD sum rules at finite density. M.P. is grateful to the McGill university nuclear theory group, where this work 
was started, for the hospitality extended to him during his visit. The work of R.H. was funded by a fellowship 
of Deutscher Akademischer Austauschdienst (DAAD). The work of A.K.D.M. has been supported by the Natural 
Sciences and Engineering Research Council of Canada and the Fonds FCAR of the Quebec Government. The work 
of M.P. was supported in part by the Department of Energy under Grant No. DE-FG02-94ER40823. 



[1] G. Q. Li, C. M. Ko and G. E. Brown, Phys. Rev. Lett. 75, 4007 (1995), H. Shiomi and T. Hatsuda, Phys. Lett B 334, 281 
(1994), H. C. Jean, J. Piekarewicz and A. G. Williams, Phys. Rev. C 49 1981 (1994), A. K. Dutt-Mazumder Nucl. Phys. 
A 611, 442 (1996) 

[2] G. E. Brown and M. Rho, Phys. Rev. Lett 66, 2720 (1991) 

[3] B. M. Preedom and G. S. Blanpied, CEBAF Report No. PR-89-001, 1989; M. Masera for the HELIOS collaboration, Nucl. 
Phys. A 590, 93c (1995), J. Friese et al, Prog. Part. Nucl. Phys. 42, 235 (1999); M. Kossov et ai, TJNAF proposal No. 
PR-94-002 (1994) 

[4] M. S. Shifman, A. I. Vainshtein, and V. I. Zakharov, Nucl.Phys.B147, 385 (1979). 

[5] M. S. Shifman, A. I. Vainshtein, and V. I. Zakharov, Nucl.Phys.B147, 448 (1979). 

[6] M. S. Shifman, A. I. Vainshtein, and V. I. Zakharov, Nucl.Phys.B147, 519 (1979). 

[7] A. I. Bochkarev and M. E. Shaposhnikov, Nucl.Phys.B268, 220 (1986). 

[8] T. Hatsuda and S. Houng Lee, Phys.Rev.C46, 34R (1992). 

[9] T. Hatsuda, S. Houng Lee, and H. Shiomi, Phys.Rev.C52, 3364 (1995). 



8 



[10] M. Asakawa and C. M. Ko, Nucl.Phys.A560, 399 (1993). 

[11] X. Jin and D. Leinweber, Phys.Rev.C52, 3344 (1995). 

[12] Y. Koike, Phys.Rev.C51, 1488 (1995). 

[13] Y. Koike, A. Hayashigaki, Prog.Theor.Phys.98, 631 (1997). 

[14] S. Leupold, W. Peters and U. Mosel, Nucl.Phys.A628, 311 (1998). 

[15] F. Klingl, N. Kaiser and W. Weise, Nucl.Phys.A624, 527 (1997). 

[16] A. K. Dutt-Mazumder and B. Dutta-Roy, Phys.Lett.B399, 196 (1997). 

[17] W. Broniowski and W. Florkowski, Phys. Lett B 440,7 (1998) 

[18] A. Bhattacharyya, Mod.Phys.Lett.A13, 2585 (1998). 

[19] A. Zhitnitsky, Phys.Rev.D55, 3006, (1997). 

[20] There are several estimates of 7 as quoted in T. Hatsuda, E. M. Henley, Th. Meissner, and G. Krein, Phys. Rev. C49, 452 
(1994). 

[21] M. Chick, E. Reya, and A. Vogt, Eur.Phys.J.C5, 461 (1998). 

[22] L.M. Barkov et al, Nucl.Phys.B256, 365 (1985). 

[23] W. Peters, M. Post, H. Lenske, S. Leupold, and U. Mosel, Nucl.Phys.A632 109 (1998). 



9 



